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STELLAR POPULATION VARIATIONS IN THE MILKY WAY'S STELLAR HALO 
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ABSTRACT 

If the stellar halos of disk galaxies are built up from the disruption of dwarf galaxies, models predict highly 
structured variations in the stellar populations within these halos. We test this prediction by studying the ratio 
of blue horizontal branch stars (BHB stars; more abundant in old, metal-poor populations) to main- sequence 
turn-off stars (MSTO stars; a feature of all populations) in the stellar halo of the Milky Way using data from 
the Sloan Digital Sky Survey. We develop and apply an improved technique to select BHB stars using ugr color 
information alone, yielding a sample of ~9000 g < 18 candidates where ~ 70% of them are BHB stars. We 
map the BHB/MSTO ratio across ~ 1/4 of the sky at the distance resolution permitted by the absolute magnitude 
distribution of MSTO stars. We find large variations of BHB/MSTO star ratio in the stellar halo. Previously 
identified, stream-like halo structures have distinctive BHB/MSTO ratios, indicating different ages/metallicities. 
Some halo features, e.g., the low-latitude structure, appear to be almost completely devoid of BHB stars, whereas 
other structures appear to be rich in BHB stars. The Sagittarius tidal stream shows an apparent variation in 
BHB/MSTO ratio along its extent, which we interpret in terms of population gradients within the progenitor 
dwarf galaxy. Our detection of coherent stellar population variations between different stellar halo substructures 
provides yet more support to cosmologically motivated models for stellar halo growth. 

Subject headings: galaxies: evolution - 
stellar content 



- galaxies: bulges — galaxies: spiral — galaxies: general — galaxies: 



1. introduction 



Over the last decade, it has become clear that a sizeable, per- 
haps dominant, part of the stellar halo of the Milky Way is com- 
posed of the accumulated debris from the disruption of distinct 
satellite galaxies. Starting with the discovery and characteri- 
zation of the ti dal debris from the disruption of the Sagittarius 
dwarf galaxy (Ibata et al. 1995; Majewski et al. 1 120031) . a vari- 
ety of other stellar ha lo overdensities have been discovered : the 
low-latitude stream (lYannv et al] 120031: llbata et al] 12003b. the 
Virgo and Hercu l es-Aquila overdensities (Duffau et al] 120061: 
lJuric et al.1 120081: iBelokurov etaD l2007ah. and a number of 
smaller tidal streams l Odenkirchen et al. 2001; Grillmair 2 0061: 



Grillmair & Dionatos 2006; Grillmair & Johnson 2006; Belokurov 
2007bl:lGrillmaid2009l:lNewberg. Yannv & Willettl2009l) . Mod- 
els of stellar halos that formed in a cosmological context through 
the di sruption of dwarf galaxies have been developed and re- 
fined [ d BullQck et a l. 2001; B ullock & Johnstod2005trAbadi et al.1 
120061: [Cooper et al] 120101) . providing quantitative predictions 
for the structure and stellar content of stellar halos formed in 
such a fashion. Remarkably, the structure and degree of sub- 
structure of the model stellar halos appear to be in quantitative 
agreement with observations of the Milky Way's stellar halo 
dBullock & Johnstonll2005llBell et alJl2008h : stellar halos with 
~ 1O 9 L , a roughly r~ 3 density profile, and showing rich sub- 
structure on a variety of spatial scales. 

A clear prediction of such a picture is that any population 
differences in the stellar halo, arising from different progeni- 
tor populations, should have a similar morphology to the stel- 
lar density inhomogeneities. The metallicity-mass correlation 
of satellite galaxies will translate into variations in stellar halo 
metallicity, as debris from larger satellites will have higher 



metallicity than those of lower-luminosity satellites or globu- 
lar clusters. Furthermore, one expects to see distinctive signa- 
tures in age and detailed elem ent abundance patterns, reflecting 
when satellites were accreted (iRobertson et al.ll2005LlFont et al] 
120061: see also lTumlinsonl2010l) . Signatures of population inho- 
mogeneity have been detected: hints of ch emical signatures o f 
kinematically-detected local streams (e.g.. iHelmi et al] 12006), 
small variations in metallicity within the nearby parts of the 
stellar halo, especially associated with the low-latitude stream 
(e.g. Ilvezic et all [200 8 fl, the tenden cy for the halo to have 
metal-poorer stars at larger radii (e.g. JCarollo et al.|[2007l) . and 
a recent color-magnitude diagram (CMD) fitting analysis of 
■^-SJipan Digital Sky Survey (SDSS) stripes showing both a radial 
metallicity gradient a nd large density and m etallicity variations 
around that gradient (Ide Jong et all l20ldb . Yet, for the most 
part, the observational signatures of such population variations 
are accessible for relatively bright stars only (high S/N ugriz 
imaging, or moderate S/N spectra). Accordingly, those radii in 
our own stellar halo where substructure is rich (10 < r < 40 kpc) 
are not probed particularly well. Some of the most impressive 
evidence for stellar halo population variations are from M31, 
where significant variations in stellar age and metallicity hav e 
been documen ted dBrown et ap l2006: Ric hardson et al Il2008|) 
In particular. llbata et al.l (l2007h and lMcConnachie et al.1 (l2009h 
have used the color of the red giant branch stars (practically 
inaccessible in the Milky Way's halo because distances are un- 
certain) to provide panoramic views of the stellar populations in 
the M31 halo, showing clearly that stellar population variation 

1 A particular concern in this case is that the low-latitude stream may have a 
significant or even dominant contribution from the debris torn off of the Milky 
Way disk; thus, the detection of a distinctive population may result from disk 
debris superimposed on a relatively homogeneous stellar halo. 
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can be associated with particular substructures. 

In this paper, we present an exploration of stellar popula- 
tion variations and their morphology within the stellar halo of 
the Milky Way. We make use of the ratio of blue horizontal 
branch (BHB) stars to main-sequence turn-off (MSTO) stars as 
our main observational diagnostic of stellar population varia- 
tions. While the exact conditions under which BHB stars de- 
velop are not completely clear (post-main sequence- stars with 
core helium burning and hydrogen s hell burning, with rela- 
tively low hydrogen envelope masses; iHoyle & Schwarzschildl 
Il955|) . studies of globular clusters demonstrate that they sign- 
post particularly ancient and metal-poor stellar populations. We 
take MSTO stars (present in all populations) as a proxy for the 
general stellar content of the halo; thus, variations in BHB to 
MSTO ratio qualitatively probe the ratio of ancient and metal 
poor stars to all stars. Despite our lack of a complete under- 
standing of which populations harbor BHB stars, these stars are 
in many other respects ideal probes of population variations in 
the stellar halo: they can be selected from the SDSS imaging 
and spectroscopic dataset, they are relatively luminous (allow- 
ing probing to large distances), and they are nearly standard 
candles. 

We outline and apply an improved color selection method 
to isolate a sample of over 9000 g < 18 high-probability BHB 
star candidates from SDSS data (^2]). The MSTO stars are 
color- selected in a simple empirical fashion to encompass the 
full range of MSTO colors in the stellar halo; we therefore use 
MSTO stars as a proxy for the general stellar content of the stel- 
lar halo. We then explore the relative distributions of BHB to 
MSTO stars in the Heliocentric distance range 5 < r/kpc < 30 
across ~ 1/4 of the celestial sphere, providing a panoramic 
view of stellar population variations in the Milky Way's stel- 
lar halo (Oh Owing to our uncertainty in the exact conditions 
under which BHB stars form, such a view of stellar population 
variations is necessarily qualitative in nature, but has the advan- 
tage that stellar population variations can be traced out to larger 
distances than is currently possible wit h complemen tary tracers 
(e.g., the lower main sequence stars of llvezic et al.ll2008l) . We 
then compare qualitatively with models of stellar halos built en- 
tirely through the disruption of dwarf galaxies (§0|), and explore 
some other implications of our results (§0). 

2. BHB AND MSTO STAR SAMPLES 

We use publicly available catalogs derived from th e Sloan 
Digital Sky Survey's (SDSS) Data Release Seven (DR7: lAbazaiian 
2009) to select BHB and MSTO stars for further study. The 
SDSS is an imaging and spectroscopic survey that has so far 
mapped ~ 1/4 of the sky. Imaging data are produced si- 
multa neously in five photometric bands, namely u, g, r, / , 
and z (iFukugita et al.ll 19961 : iGunn et al.lfl998l: iHogg et al.ll200ll: 
iGunn et al.ll2006l) . The data are processed throug h pipelines to 
measure photometric and astrometric properties dLupton et al] 
1999HStoughton et al]l200ilSmith et aD l2002l:[Pier et all2003t 
Ivezic et al.ll2004l; iTucker et a l. 2006) and to select targets for 
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Fig. 1. — Illustration of the recovery of BHB stars for the globular cluster 
NGC 5024. The bottom panel shows the cluster CMD over a wide range of 
colors and magnitudes; the top panel zooms into the area richest in BHB stars. 
Stars selected as BHB candidates by -0.5 < g-r < and 0.8 < u-g < 1.6 
are shown as crosses, while those with > 50% probability of being a BHB star 
as estimated by our method are marked with an additional diamond. 



refer to these as a blue plume, e.g., iBellazzini et al] [2006a). 
Fig. Q] illustrates this issu^E showing BHB stars at g ~ 16.85 



spectroscopic follow-up dBlanton et al.l2003l;IStrauss et al. 12 002). 

To obtain full area coverage and maximal distance range, we 
must rely on BHB star samples selected by ugr color infor- 
mation alone. One of the well-known practical challenges of 
using BHB stars for stellar halo science is that genuine BHB 
stars are outnumbered by in trinsically fainter blue straggler 
stars (as discussed at length b y Kinman, Suntzeff & Kraft ] 1994 



Wilhelm, Beers & Grav 1999, Clewlev et al. 2002, ISirko et al] 
2004 Kinman et al. 2007, and Xue et al. 2008; some authors 



et gnd less luminous blue straggler stars with similar colors at 
T8 < g < 20. These blue straggler stars are a common fea- 
ture of globular clusters, dwarf galaxies and the stellar halo 
(e.g. ISandagdll953l: IBellazzini etaDl2006al: fPreston & Snedenl 
l2000l: ICarney et al .11200 ll) . Blue stragglers are somewhat enig- 
matic, as they are intermediate mass main sequence stars in 
(old) populations where the MSTO is at lower masses. It is 
thought that there are at least two po ssible formation r outes: 
stellar collisions in dense star clusters (Hills & Dav 197 q), and 
mass transfer in the evolution of cl ose binaries (McCrea 11964 
or potentially in triple star systems : |Perets & Fabry ckyl2009l) in 
low-density environments (where the latter low stellar density 
route is more relevant here; see, e.g.. JPreston & Snedenl [20001 
iMathieu & Gelledl2009llKnigge. Leigh & Sillsll2009l) ~ 

2 Because NGC 5024 is particularly BHB star rich, the blue stragglers do 
not dominate in this figure. When translated to the global halo population, 
and accounting for the stellar halo density profile and different distances of 
BHB and blue straggler stars at a given magnitude, BHB stars are typically 
outnumbered by a significant factor by blue stragglers. 
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Fig. 2. — Illustration of the BHB star selection method which identifies 
the optimal BHB selection region in the u-g vs. g-r color plane using exten- 
sive spectroscopic identification using the SDSS. Only stars with spectroscopic 
identifications are shown. Stars that are argued to be BHB stars on the basis of 
their spectra by Xue et al. ( 2008) are shown in blue. Contours show the frac- 
tion of candidates determined to be BHB stars (contours are 10%, 30% and 
50%; the 50% contour is thicker than the others). In what follows we consider 
BHB stars within the 50% contour only. Overplotted are our appro ximation to 
the BH B color selection limits (red), and the 'stringent' color cut of Sir ko et all 
(12001 green). 



2.1. Color selection of BHB stars 

Choosing stars in the broad color cut broad cuts -0.5 < 
g-r < and 0.8 < u-g < 1.6 gives a sample that is >3:1 
blue stragglers t o BHB stars for the full g < 19.5 sample of 
IXue et al.1 (120081) . and 60% blue stragglers for the g < 18 sub- 
sample of interest in this paper. This ratio is based entirely 
on the sample of stars with SDSS spectroscopy as reported in 
IXue et al] (12008b . The algorithm used to select stars for SDSS 
spectroscopy changed a number of times during the survey, and 
does not select uniformly in color-color space; these contam- 
ination fractions are affected by this selection (we present al- 
ternative estimates of contamination fractions derived in other 
ways later). Yet, it is clear that such a level of contamination is 
a significant challenge for quantitative analysis of the BHB star 
population using photometry alone. 

Our goal here is to combine the extensive SDSS spectroscopy 
with SDSS photometry to devise color selection criteria that are 
significantly improved in their trade-off between BHB purity 
and sample size (similar in spirit to the 'strin gen t' color cut of 
Sirko et al.l2004t see also lClewlev et al .120051 and lKinman et al.l 
20071 ? or selection of BHB stars from ugr photometry). Us- 
ing a calibration sa mple of BHB star candidates with spec- 
tra (IXue et al.l l2008. following the methodology of Sir ko et all 
12001 7 we have determined the probability that a BHB candi- 
date with given u-g and g-r color is a spectroscopic BHB 
star; the rest of the candidates are other blue stars, and are pre- 
sumably mostly blue stragglers. Fig. [2] illustrates this process, 
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Fig. 3. — BHB star candidates in the direction of RA/Dec=(204,5). The 
lower panel shows the logarithmic density of stars in the CMD in gray scale 

o 

in an area of radius 4 , with those stars identified as BHB candidates on the 
basis of their ugr colors being shown with crosses. The upper panel shows 
the distance modulus distribution of BHB candidates along this line of sight, 
showing two overdensities; the more distant of these at m-M ~ 18.3 is part of 
the tidal tail from Sagittarius. It is worth noting that blue straggler contamina- 
tion has a much broader absolute magnitude distribution; the confinement of 
this overdensity to such a small apparent magnitude range is further evidence 
that the sample is dominated by BHB stars (that have a small range in absolute 
magnitude). 



showing stars with g < 18 (appropriate for BHB stars in the He- 
liocentric distance range considered here) and with spectra from 
the SDSS ( Xue et al]|2008l) . Blue data points show stars that are 
very likely to be BHB stars on the basis of their spec tra (a con- 
tamination of much less than 10% has been argued by lXue et~ai] 
120081 and ISirkoetaD 120041 at these g < 18 limits). The thick 
contour outlines the region of color-color space where the frac- 
tion of BHB candidates that are spectroscopically-classified 
BHB stars is > 50% and there were more than 16 stars in a bin 
of 0.025 x 0.04 mag (to ensure a high enough number of stars to 
measure the BHB star fraction with some fidelity). This region 
is well-approximated by the selection region shown in red in 
Fig. [2 0.98 < u-g < 1.28, -0.2 < g-r < -0.06 and excluding 
the region with ([w-g-0.98]/0.215) 2 +([g-r+0.06]/0.17) 2 < 
1 (for mulated similarly to the 'stringent' color cut of lSirko et al] 
l2004l shown in green). 

2.2. Completeness and contamination of color-selected 
samples 
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We first quantify the likely completeness and contamination 
of the color- selected g < 18 sample that is the main focus of 
this paper. We estimate the completeness and contamination 
in two ways: using the spectroscopic sample, and independent 
spot checking on BHB star-rich globular clusters. 

Our primary estimate of contamination and completeness 
comes from the spectroscopic training sample used to deter- 
mine the color cuts. Given the average spectroscopic complete- 
ness within the > 50% contour, and weighted by the number 
of stars in our photometric catalog per bin in color-color space 
(this to a certain extent corrects for the non-uniform color se- 
lection of stars with SDSS spectra), we expect ~ 67 - 74% of 
the BHB candidates within the selection region to be real BHB 
stars (where the two ends of the range assume that 90% to 100% 
of the spectroscopic BH B stars are act ually BHB stars, follow- 
ing lSirko et al.l2004l and lXue et al.ll2QQ8h . Furthermore, we find 
that the color cuts select ~ 60% of the BHB stars in the spec- 
troscopic sample. 

Secondary estimates of contamination and completeness are 
possible through the investigation of BHB stars in globular 
clusters and dwarf galaxies. This method has the advantage 
that an independent indicator of the nature of the stars is avail- 
able through their position on the CMD, but have the disad- 
vantage that there are very few clusters in the relevant distance 
range, and that the number of BHB stars per cluster is typ- 
ically very modest (dwarf galaxies were also tested, but are 
more distant, and give results qualitatively consistent with the 
picture presented here but are of limited quantitative use). In 
Fig. [TJ we illustrate this method to the BHB star-rich globular 
cluster NGC 5024. Stars selected as BHB candidates by the 
broad cuts -0.5 < g-r < and 0.8 < u-g < 1.6 are shown as 
crosses, while those with > 50% probability as estimated by our 
method are shown as diamonds. There are 30 actual BHB stars 
in this sample at the distance of NGC 5024 (defined as having 
-0.4 < g-r < and 16.6 < g < 17.1); our method recovers 
11 of them with only 2 contaminants (i.e., a completeness of 
~ 37% and a contamination of ~ 15%, consistent with our ex- 
pectation of < 30% contamination). There are indications that 
NGC 5024 may have an unusually low fraction of recovered 
stars; NGC 5053 has 5/7 (~ 70%) BHB stars recovered (with 
no contaminants). 

Considering both tests together, we conclude that this proce- 
dure, when applied to g < 18 stars with SDSS-quality photom- 
etry, gives a sample of ~ 50% of the real BHB star population 
with < 30% contamination, a considerable improvement over 
the ~60% contamination given by much more simplistic pho- 
tometric cuts. 

While in this paper we restrict our attention to BHB stars 
with g < 18, we characterize the likely performance of these 
cuts for fainter samples with colors characteristic of BHB stars. 
To estimate the run of photometric uncertainty as a function of 
measured apparent magnitude in u, g and r-bands , we used re- 
peat imaging of the SDSS stripe 82 as tabulated by lBramich et al.l 
( 2008). We then imposed the typical uncertainties of stars with 
BHB colors at g ~ 19 and g ~ 20 on the sample of real g < 18 
stars used to derive the BHB star selection, and checked how the 
fraction of the real BHB stars that are recovered, and the con- 
tamination by BS stars, should depend on apparent magnitude. 
Recall for the g < 18 stars, the recovery/contamination fractions 
for the spectroscopic sample are (60%, 25%), i.e., the cut iso- 
lates 60% of the population with only 25% contamination. For 
g=l9 and g = 20, our estimates of the recovery /contamination 



fractions are (40%, 40%) and (20%, 50%) respectively, i.e., at 
g = 19 the situation is improved over the 60% contamination 
in the broad -0.5 < g-r < 0, 0.8 < u-g < 1.6 color bin by 
applying the cuts described here, but at g = 20 the cuts throw 
out most of the genuine BHB stars without improving contam- 
ination significantly. These fractions are consistent with our 
investigations of the performance of cuts on more distant BHB- 
rich targets, e.g., the Sextans or Draco dwarf galaxies (1/4 of 
the BHB stars recovered, at g ~ 20). 

2.3. Estimates of BHB star absolute magnitude 

We estimate the absolute magnitudes for the BHB star candi- 
dates by in terpolating the g- r-dependent M g calibration from 
Table 2 of fsTrko et alJ (120041) with [Fe/H]= -1 (the [Fe/H]= -2 
iSirko et al. 1 120041 calibration is different by < 0.03 mag). Tests 
of the quality of this calibration are encouraging; the rms of 
the color-derived absolute magnitudes of BHB stars in globu- 
lar clusters is ~ 0.1 mag (a measure of the absolute magnitude 
scatter at a single metallicity), and the overall distance modulus 
estimates of the clusters and dwarf galaxies agree with those 
determined by independent means to < 0.1 mag (a rough con- 
straint on the degree of metallicity dependence of the BHB ab- 
so lute magnitude calib ration; these results are reported in full 
in lRuhland etaDl2010l) . 

We also compared the cali bration of BHB absolute mag- 
nitude from ISirko et all ( 20041) with the m etallicity -dependent 
BHB star models from iDotter et al.1 (120071) and [Potter et al.l 
(2008). These models m ake use of more re alistic input physics 
than the modeling from Sir ko et al.l (|2004|); however, the ugr 
colors of BHB stars from lDotter et al.l (120081) do not match the 
colors of SDSS BHB stars particularly well (they are offset in 
particular in u-g color by ~ 0.1 mag, and may be offset by a 
smaller amo unt in g - r color; in contrast, the more rudimentary 
modeling by lSirko et al.ll2004l was designed to match the colors 
well). These offsets in color prohibit a quantitative assessment 
of which stars would fall in the selection region; selecting in- 
stead those BHB stars reddest in u- g at a given g- r (a qualita- 
tive version of our selection method), corresponding to metal- 
licities of [Fe/H] < -1.5, gives a full range in AM g < 0.3 mag 
at a given g-r. Taken together, the tests of absolute magnitude 
cal ibration on c l usters and dwarf galaxies in conjunction with 
the lDotter et al.l (120081) metallicity-dependent BHB star models 
paints a picture in which the absolute magnitudes adopted for 
our BHB candidates are accurate to < 0. 1 in overall zero point, 
and will scatter < 0.2 mag at worst - far less than the ~ 0.9 mag 
absolute magnitude scatter of our MSTO sample. 

The final sample of high-probability BHB stars selected us- 
ing the 50% color contour in Fig. [2 and with g < 21 is more 
than 14000 stars in the North Galactic Cap selection region we 
use in this paper (b > 20 for 180 < / < 240, and b > 30 for 
/ < 180 and / > 240, dictated largely by the SDSS coverage), 
and more than 9000 stars in the Heliocentric distance range 5- 
30kpc (corresponding approximately to g < 18). 

2.4. Selection of MSTO stars, and comparison of MSTO and 

BHB stars 



Following lBell et al.l (120081) . we select over 4.3 million MSTO 
stars in the foreground extinction-corrected color range 0.2 < 
g-r < 0.4 and 18 < r < 22; this color and magnitude range was 
selected empirically to encompass the most densely-populated 
bins of color s pace for the halo MSTO stars (see Fig. 1 of 
Bell et al.ll2008h . Examination of globular clusters and model 
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CMDs demonstrates that MSTO stars selected in this fashion 
have a roughly Gaussian distribution of absolute magnitudes 
wi th a mean M r ~ 4.5 and a rms cr magi MSTO = 0-9 (see Fig. 2 
of iBell et"ai1l2008l) . By choosing a particular color range for 
MSTO star selection, we are in principle susceptible to changes 
in the main- sequence stellar population from location to loca- 
tion in the stellar halo (through variation in the metallicity and 
age). In practice, this turns out to be a minor concern, as the 
color of the stellar halo MSTO (the blue edge of the turn-off in 
particular) varies little across the SDSS area; the stellar popu- 
lation in the stellar halo varies little enough that it is essentially 
undetectable using MSTO color alone. 

Fig. [3] illustrates the properties of the stellar sample in a 
particular line of sight towards the Sagittarius tidal tail. This 
shows the CMD of (primarily halo) stars in a 4° radius cone 
around (RA,Dec)=(204,5). MSTO stars in the halo are very nu- 
merous, having 0.2 < g-r < 0.4 and g > 18 (brighter MSTO 
stars with g < 18 have g-r ~ 0.4, representing more metal- 
rich thick disk stars). BHB candidates are selected by virtue 
of their ugr colors using the method illustrated in Fig. O and 
are shown as crosses. Because BHB stars have a small scat- 
ter in absolute magnitude, BHB star-rich halo overdensities (at 
m-M ~ 16.5 and m-M~18.3)are clearly-defined. The corre- 
sponding MSTO overdensities, about 4 magnitudes fainter, are 
not anywhere near as clearly defined owing to the larger spread 
in MSTO absolute magnitudes. In this paper, we will be analyz- 
ing the two populations in concert as an easily-measurable but 
largely qualitative diagnostic for stellar population variations in 
the stellar halo of the Milky Way. 

3. POPULATION VARIATIONS IN THE STELLAR HALO OF THE 
MILKY WAY 

Analyzing the ratio of BHB to MSTO stars requires that we 
compare the stars in the same 'effective volume'. We con- 
volve the BHB star distance moduli with a Gaussian with width 
o 'mag, msto = 0.9 mag (referred to hereafter as the degraded BHB 
sample) to yield the same scatter between degraded distance 
modulus (m-M) degraded an d true (jn-M)tme as is estimated for 
the MSTO stars. Then we can compare the 'distance smoothed' 
distribution of BHB stars, with (m- M) degraded , directly and 
quantitatively to the distribution of MSTO stars with the same 
(m-M) M sTO range (i.e., BHB stars are compared with MSTO 
stars that are ~4 mag fainterfl The results are illustrated in Fig. 

13 

Fig. [4] displays various aspects of the comparison between 
BHB and MSTO star distributions, and we work through them 
in turn. The center column of panels shows the distribu- 
tion of MSTO stars in four distance modulus slices: 13.5 < 
m-M < 14.5, 14.5 < m-M < 15.5, 15.5 < m-M < 16.5, and 
16.5 < m-M < 17.5, or distances of 5-8, 8-13, 13-20 and 
20 - 32 kpc respectively, smoothed in angular position with a 
Gaussian kernel with a = 0.5°. The gray scale shows the num- 
ber of stars in 0.5° x0.5° bins with a linear mapping between 
intensity and stellar density. Black areas lack SDSS DR7 cov- 
erage. These panels show various pre viously described fea- 
tures of the Milky Way stellar halo (e. g . , iBelokurov et al.l 120061: 
iNewberg et ai1l2007l: Bell et alJ l2008): superimposed on a rel- 
atively smooth distribution of stars is the Sagittarius stream 

3 It should be noted that this shift in magnitude is important also when con- 
sidering blue straggler contamination. Blue stragglers are 1-3 mag fainter than 
BHB stars, therefore any blue straggler star contamination should show up 
at larger equivalent distance modulus (because the blue stragglers are being 
treated as brighter BHB stars). 



fe^ llbataet1u r lll995l: iMaiewski etaD 120031: IBelokurov etaD 
2006) stretching from the Galactic anticenter (m-M ~ 16) to 
the Galactic center (m-M ~ 17). Near the North Galactic Pole 
(NGP) at m -MS, 16 is the diffuse V i rgo overdensity (e . g . , 
iDuffau et alJ l2006kl iNewberg et al.ll2007l: lJuric et al.ll2008l) . At 
m-M < 15 in the Galactic anticenter direction at low galac- 
tic latitudes is the low-la titude structure (also referred to as the 
Monoceros stream; e.g.jNewberg et aL 2002; Pena rrubia et al.l 
120051: iMomanv et al.ll2006l) . At / ~45 ° and b ~ 45° is t he dif- 
fuse Hercules-Aquila structure (IBelokurov et al.ll20 07a). One 
can also see a variety of less prominent features: globular 
clusters as 'hot pixels', and at m-M ~ 17 stretching from 
(Lb} ~ (250,50) to (l b) ~ (160 40) i s the Orphan Stream 
( Belokurov et ai1l2007bl: lGrillmaidl2006l) . 

In the left-hand column of panels in Fig. 01 we show the BHB 
star distribution in the same distance bins, and convolved with 
a distance modulus kernel of width ctmsto = 0.9 mag to match 
the MSTO absolute magnitude scatter; an angular smoothing 
of a = 1° is applied to suppress the much larger shot noise 
due to the 100-fold smaller sample. Thus, the left-hand pan- 
els reflect BHB star maps of the same halo distance slices 
as the MSTO panels. One can see some similarities, as well 
as some striking differences, between the distributions of the 
BHB and MSTO stars. These differences are quantified and 
illustrated in the right-hand panels, which show a color rep- 
resentation of the BHB to MSTO star ratio. Here, the distri- 
butions of BHB and MSTO stars are smoothed by a Gaussian 
kernel with an angle of a = 6° to increase S/N. A bluish-white 
color denotes BHB/MSTO ratio of > 1/50 and red denotes 
BHB/MSTO < 1/100. The intensity scales with the number 
of stars (dark shades show regions with few stars, whereas light 
shades denote regions with many stars). Poisson noise is signif- 
icant on small scales, but contributes less than 3% on the scale 
of large structures such as the low-latitude stream. The pri- 
mary uncertainties are systematic in origin: changing the prob- 
ability threshold for a BHB star to be included in the sample 
from 50% to 30% tests the influence of BHB candidate purity 
and blue straggler contamination on the results. Such a change 
produces differences of ~ 10% in relative BHB/MSTO ratio 
(i.e., variations in the BHB/MSTO from place to place on the 
mapfl Another source of systematic error is possible varia- 
tions in MSTO ab solute mag nitude or scatter as a function of 
position. Fig. 2 of lBell et alJ d2008h shows that there are some 
variations in MSTO absolute magnitude and rms from globular 
cluster to globular cluster, driven by differences in stellar pop- 
ulation age/metallicity. Such variations may be present in the 
stellar halo. We have tested the influence of potential variations 
in MSTO absolute magnitude (a < 15% effect in BHB/MSTO 
ratio for variations in MSTO absolute magnitude of 0.5 mag or 
less) and scatter in the luminosity of MSTO stars (a few percent 
effect in BHB/MSTO ratio for 20% changes in the assumed 
MSTO scatter). 

The right set of panels in Fig. [4] demonstrate clearly the 
key result of this paper: there are significant variations in the 
BHB/MSTO ratios in the stellar halo of the Milky Way. Per- 
haps more importantly, these variations coincide spatially with 
different known halo structures, as expected if the various over- 
densities were composed of the debris from the disruption of 

4 Of course, the absolute BHB/MSTO ratio changes by 30%, as the number 
of BHB star candidates included in the analysis is increased by reducing the 
BHB probability threshold; such a change in absolute ratio is not important for 
our purposes. 
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individual progenitors. Some prominent structures (e.g., the 
low-latitude structure) appear to be almost completely devoid 
of BHB stars (as illustrated by the red band at low galactic lati- 
tudes in the Galactic anticenter direction in the upper right-hand 
panel of Fig. H BHB/MSTO - 1/90). Other structures, e.g., 
the Virgo overdensity at / ~ 300° and b ~ 70° or the Hercules- 
Aquila overdensity at / ~ 45° and b ~ 45°, appear to be consid- 
erably richer in BHB stars (both structures have BHB/MSTO 
~ 1 /40). Such differences are highly significant; recalling that 
the primary uncertainties are systematic, we find that Virgo and 
Hercules- Aquila are a factor of 2 to 2.5 richer in BHB stars 
than the low-latitude structure (where the range includes the 
effects of the random and systematic uncertainties discussed 
above). Strikingly, parts of the Sagittarius tidal stream (near 
the NGP and at m-M ~ 17) appear to be relatively rich in 
BHB stars (BHB/MSTO - 1 /55), while the part at m-M ~ 16 
and / ~ 200° and b ~ 50° is considerably poorer in BHB stars 
(BHB/MSTO ~ 1/80), implying a population gradient along 
the stream. 

Fig- S] also shows that contamination of the BHB sample by 
blue stragglers would not affect our key result. Areas with al- 
ready low BHB/MSTO ratios (the low-latitude stream and parts 
of Sagittarius) would imply even lower BHB/MSTO in the pres- 
ence of blue straggler contamination. Furthermore, even for 
structures such as Virgo or Hercules-Aquila, the association 
of BHB stars to particular structures already present in MSTO 
stars (and their having characteristic BHB/MSTO ratios) im- 
plies that blue straggler contamination has a small effect, at 
least for m-M < 17.5. The only possible exception to this 
is the BHB star-rich section of Sagittarius, which lies in a sim- 
ilar area of the sky and 1-2 magnitudes more distant than the 
BHB star-rich Virgo overdensity. It is possible that some of 
these BHB stars are blue stragglers from Virgo, but there are 
two arguments that blue straggler contamination is not domi- 
nant. Firstly, the morphology of the BHB stars in the ~ 27 kpc 
slice in Fig. 01 argues strongly that most of the BHB star can- 
didates are in fact BHB stars associated with Sagittarius — the 
distribution of BHB star candidates is elongated like the Sagit- 
tarius tidal stream, and extends past / = 270 where there are no 
stars in the Virgo overdensity (as can be seen in the ~ 1 1 kpc 
slice). Secondly, there are numerous blue stars along those lines 
of sight that are 1-2 mag fainter than the BHB stars we identify 
here; we interpret these numerous fainter blue stars as the blue 
straggler content of Virgo an d S agittarius (this is ill ustrated in 
iNiederste-Ostholt et alJl201Qi and lRuhland et alJl2010h , 

4. COMPARISON WITH MODELS OF STELLAR HALO FORMATION 
IN A COSMOLOGICAL CONTEXT 

One motivation for this analysis was to test further the hy- 
pothesis that the stellar halo is composed of the debris from the 
disruption of dwarf galaxies. Cosmologically- motivated simu- 
lations of stellar halo formatio n such as those of Bull ock & Johnston! 
(2005]) or lCooper et al] (120101) allow reaso nably direct compar- 
ison o f the observations and simulations. Bull ock & Johnston] 
(120051) have publicly releasecQ the results of their simulations, 
including the luminosities, ages, metallicities, and element abun- 
dance s of star particles (see also lRobertson et alj2005tlFont et al] 
2006) in eleven simulated stellar halos. Our main observational 
diagnostic, the spatially-varying ratio of MSTO to BHB stars, 
cannot be predicted directly by the models: while old and metal 
poor populations contain a significant population of BHB stars, 

5 http : / / www .astro . Columbia . edu/^kvj / ha los /| 



it is not completely clear what the exact conditions are under 
which stellar populations are rich in BHB stars. An illustra- 
tion of this challenge is the second-parameter problem: glob- 
ular clusters with very similar metallicities can have very dif- 
ferent horizontal branch morphologies, where it is unclear what 
parameter or parameters drive these va riations (e.g., age or He- 
lium content; see, e.g., ICatelanl [2009 for a discussion of this 
issue). Obviously, this uncertainty makes a correct BHB pre- 
scription impossible to implement. 

Bearing in mind the practical difficulty of predicting the ra- 
tio between BHB and MSTO stars, we probe instead a metric 
that is better-defined in the models that gives some qualitative 
idea of population variations in stellar halos: stellar metallic- 
ity. We present in Fig. [5] a map of r-band luminosity-weighted 
stellar metallicity of three si mulated stellar h alos (Model num- 
bers 4, 2, and 1 of Fig. 13 of iBell et al.ll2008l corresponding to 
halo n umbers 8, 5 and 2 respectively from lBullock & Johnston! 
2005) in the same distance bins as the observations (and con- 
volved with the same absolute magnitude scatter of 0.9 mag as 
we use for the observations in Fig. [4]). In order to produce 
this map, we have smoothed the most luminous star particles 
to reduce noise (although some noise in the maps is visiblefl 
The maps are luminosity- weighted (intensity scale), and color- 
coded by stellar metallicity (blue shows [Fe/H] = -1.5, red 
shows [Fe/H] = -0.5). The three halos were chosen to have 
differing amounts of substructure; model 4 (left) has a modest 
amount of substructure and model 1 (right) has prominent sub- 
structure, especially in the most distant bins. The key point to 
take away from Fig. [5] is that it clearly demonstrates the predic- 
tion that different substructures should have distinctive stellar 
populations (in this case parameterized by metallicity). This is 
in qualitative agreement with the observations — a quantitative 
match will be impossible to produce until the conditions under 
which BHB stars form are substantially better understood. 

5. DISCUSSION 

We have presented a map of spatial stellar population vari- 
ations in the Milky Way's halo, using the ratio of BHB to 
MST O stars. This map is reminiscent of analogous maps for 
M3 1 dMcConnachie et all2009l) and shows that descriptions of 
stellar halo populations in terms of radial gradients alone is an 
incomplete description of the data; much of the population vari- 
ation is stream-like in morphology. 

Our results may bear on the debate about the controversial 
nature and origin of the low-latitude structure: it may represent 
debri s from a disrupted dwarf galaxy (e.g., iPenarrubia et al] 
2005), or primarily consist of m aterial stirred up from the outer 
disk of the Milky Way (see, e.g . Jlbata et al.l20 05: Mom any et al] 
2006; Kazantzidi s et al] [20081 for discussion of this possibil- 
ity). Fig. 01 shows the low-latitude structure to be almost de- 
void of BHB stars, implying in a broad sense a lack of metal- 
poor old populations. This is in contrast to many of the other 
halo substructures (at least parts of Sagittarius, Virgo, and Her- 
cules/Aquila). If it is found in the future that the outer thin disk 
of the Milky Way is deficient in BHB stars, this may provide 

6 In detail, we split star particles more lumi nous than 8 sola r luminosities 
(sufficient to produce more than 1 MSTO star; IBell et al.l f2008) into L r /8L 
particles, where L r is the luminosity of the star particle of interest. In or- 
der to determine the placement of the new particles in 3D space, we find the 
nearest particles (between 100 and 1000 particles, depending on luminosity) 
drawn from the same satellite as the particle of interest (in order to preserve 
stream/satellite morphology) and distribute the new star particles in a space de- 
fined by these nearest neighbors. In practice, the details of this procedure are 
not critical for this application. 
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tentative support to the notion that the low-latitude structure is 
composed primarily of material stirred up from off the Milky 
Way's outer disk. At the very least, the lack of BHB stars (one 
of the few options for precise distance determination) is a sig- 
nificant practical challenge for those attempting to understand 
the 3 -dimensional distribution of the low-latitude structure. 

A feature of particular interest in Fig.|4]is the change in BHB 
content al ong the Sagittarius tida l strea m. It is interesting to 
note that iNiederste-Qstholt et al.l (120101) argue for a constant 
BHB/MSTO ratio along the Sgr stream in the very region that 
we claim a signific ant difference in BHB/MSTO ratio (from 
- 1/55 to - 1/80). INiederste-Qstholt et al.l (l201Qh subtract off 
the CMD of the region with the same b but / con troi = 180- /s gr , 
then count the number of MSTO stars and BHB stars left in the 
residual CMD. We have confirmed that the BHB/MSTO ratios 
of the control areas systematically change, from ~ 1 /50 in the 
mirror region of the BHB star-rich part of Sgr, to < 1 / 250 in the 
mirror region of the BHB star-poor part of Sgr (this is apparent 
in Fig.Hl the top part of the last set of panels is relatively BHB- 
rich, whereas the lower parts of the last set of panels is poor 
in BHB stars both in Sgr and elsewhere). This change in the 
properties of the control sample drives the apparent constancy 
of the BHB/MSTO ratio in their work. It is not obvious (to us 
at least) how to properly interpret such a situation. From the 
perspective of an underlying smooth stellar halo (with abrupt 
changes in BHB/MSTO in the 'smooth' halo) with a superim- 
posed stream, subtracting off the control fields is defensible. 
Yet, from the perspective of viewing the stellar halo as a combi- 
nation of a number of structures in various stages of disruption, 
subtraction of a control field is less defensible, and would lead 
to artificial changes in the inferred properties of the structure of 
interest. 

Viewed from the latter perspective, such a difference in the 
BHB/MSTO ratio would demonstrate that population gradients 
within a progenitor galaxy could in practice lead to popula- 
tion differences in the resulting tidal debris (a gradient in the 
Sagittarius tidal tail has been d etected before using other stel- 
lar population diagno s tics; e.g.jMartmez-Delgado et al.ll2004t 
BellazzinietaD|2003 IChou et al.ll2007h . It may well be that 
the apparent abruptness of the change in the BHB/MSTO ratio 
between the Sagittarius debris in the anticenter direction and 
the NGP direction is because the debris streams in those direc- 
tions wer e stripped from Sag ittarius during different passages 
(see, e.g jLaw et al .1120051 and lLaw & Majewskill2010l for mod- 
els in which these two parts of the debris stream were stripped 
at significantly different times). Material closer to the edges 
of a satellite is preferentially stripped first; thus, in this inter- 
pretation we would hypothesize that the outermost parts of the 
Sagittarius dwarf galaxy were devoid of BHB stars, whereas the 
central parts of the dwarf were richer in BHB star^J Spectro- 
scopic follow-up of stream members, with the goal of weed- 
ing out phase-mixed stars from stars in a dynamically-cold 
stream, will help to differentiate between the two possible in- 
terpretations of the apparent p opulation gradient in Sgr (see 
iKeller. Yong & Da Costal I20T0I for such an investigation of the 
Sgr trailing arm). 

7 It is worth noting that this picture may run counter to the scenario sug- 
gested by the observed differences in e.g., BHB to red clu mp star ratio, or 
metallicity, between the Sagittarius core and the tidal streams ( Bellazzini et al. 
2006; Chou et al. 2007). In their picture the core of Sagittarius is more metal- 
rich, and the tidal tail more metal poor. Here, we are sensitive to differences 
between different parts of the tidal tail; how that relates to earlier core/tail com- 
parisons is not entirely clear. 



6. CONCLUSIONS 

In this paper, we have studied the spatial structure of stellar 
population variations in the stellar halo of the Milky Way. We 
made use of a new color selection method to isolate a sample of 
high-probability BHB star candidates, and compare their spatial 
density to that of color- selected MSTO stars (taken to represent 
the general stellar content of the halo). The abundance of BHB 
stars (vs. MSTO stars) is known to vary strongly among stellar 
populations with globular cluster-like ranges in age and metal 
abundance, where broadly speaking high BHB abundance sign- 
posts particularly old and metal-poor populations, and redder 
Horizontal Branch populations are characteristic of younger, 
more metal-rich populations, though the physical origin of such 
variations is currently debated. We mapped the relative distri- 
butions of BHB and MSTO stars across the Heliocentric dis- 
tance range 5 < r/kpc < 30 for ~ 1/4 of the celestial sphere, 
providing a panoramic view of the content of the stellar halo. 

We found large variations of the BHB/MSTO star ratio in 
the stellar halo. Most importantly, variations trace different 
previously-identified structures, indicating distinct populations 
and h ence origins for them (in common, for example, with 
M3 1 : iMcConnachie et ai1l2009l) . Some halo features, e.g., the 
low-latitude structure, appear to be almost completely devoid of 
BHB stars, whereas other structures appear to be rich in BHB 
stars. The Sagittarius tidal stream shows an apparent variation 
in the BHB/MSTO ratio along its extent, which we interpret 
in terms of population gradients within the progenitor dwarf 
galaxy leaving observable signa tures in our stellar halo. In a 
previous paper (iBell et al.ll2008l) . we had shown that the level 
of density substructure in the Milky Way's stellar halo is consis- 
tent with models (e.g jBullock & Johnstonll2005l) in which this 
component is built up exclusively from disrupted satellites. In 
this paper we have shown that another prediction of such mod- 
els is qualitatively borne out: significant population variations, 
traced by the BHB/MSTO star ratio in the Milky Way's stellar 
halo, and with a spatial structure that correlates with the den- 
sity substructures. This lends further observational support to 
the view that the stellar halo is predominantly assembled from 
the disrupted debris of dwarf galaxies. 
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Fig. 4. — Left: Map of the BHB stars in thick distance modulus slices, where the BHB distances have been 'degraded' to gm - 0.9 to match the distance resolution 
of MSTO stars. Middle: Map of MSTO stars in the same distance modulus slice, assuming a MSTO M r = 4.5. Right: Color representation of the BHB/MSTO ratio, 
smoothed with a 6 degree Gaussian. Panels are arranged in order of increasing distance modulus from top to bottom, 13.5 < m-M < 14.5, 14.5 < m-M < 15.5, 
15.5 < m-M < 16.5, and 16.5 < m-M < 17.5, or distances of 5-8, 8-13, 13-20 and 20-32 kpc respectively. In all panels, a gray scale that varies linearly with 
the star number is used. 
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Fig. 5. — A color representation of the relationship between structure and metallicity in simulated stellar halos of Bullock & Johnston ( 2005). The three columns 
show stellar halos with relatively low, average, and large amounts of substructure. The 4 rows show the simulated halos convolved with the MSTO absolute 
magnitude distribution and split into the same distance bins as used for the SDSS data: 13.5 < m-M < 14.5, 14.5 < m-M < 15.5, 15.5 < m-M < 16.5, and 
16.5 < m-M < 17.5, or distances of 5-8, 8-13, 13-20 and 20-32kpc respectively. The intensity scale is linear and scales with the number of stars; the color 
scale shows metallicity, where blue denotes [Fe/H] = -1.5 and red denotes [Fe/H] = -0.5. 



